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Abstract 

The Gamow- Teller resonance in 208 Pb is discussed in the context 
of a self-consistent RPA, based on the relativistic mean field theory. 
We inquire on the possibility of substituting the phenomenological 
Landau-Migdal force by a microscopic nucleon-nucleon interaction, 
generated from the rho-nucleon tensor coupling. The effect of this 
coupling turns out to be very small when the short range correlations 
are not taken into account, but too large when these correlations are 
simulated by the simple extraction of the contact terms from the re- 
sulting nucleon-nucleon interaction. 
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The quantum hadrodynamics (QHD) aims to describe the nuclear many- 
body system in terms of nucleons and mesons 0. Proposed initially as a 
full-fledged renormalizable quantum field theory, nowadays it is seen as an 
effective field theory, derivable, in principle, from the quantum chromody- 
namics [@. 

The relativistic mean field theory (RMFT), which can be thought as a 
mean field (Hartree) approximation to the QHD, has been applied with great 
success during the last few decades. For instance, it accounts for both i) the 
nuclear matter saturation, and ii) the ground state properties of finite nuclei 
along the whole periodic table ||. More recently, the RMFT has also been 
exploited for the description of unstable nuclei all up to the nucleon drip 
lines ||. 

Through a relativistic version of the random phase approximation (RRPA), 
various excited states and resonances have been studied in the context of the 
RMFT ||— 1| as well. Quite recently we have also reported [JTC] the first cal- 



culation of this type for the Gamow- Teller (GT) and isobaric analogue (IA) 
resonances, excited from the ground states of 48 Ca, 90 Zr and 208 Pb nuclei. 

Because of its pseudoscalar nature, the pion does not participate in the 
description of the the ground states in the RMFT. Thus, besides the nucleon 
and the Coulomb fields, only the a, uj and p mesons are usually involved in 
the calculations. Yet, in dealing with isovector excitations it is essential to 
include, together with the p meson, the 7r meson as well. This has already 
been done in our previous work flO|| , with the pseudovector pion- nucleon 



coupling / w fixed at its experimental value. For the remaining mesons, only 
the nonderivative couplings to the nucleon were included, as usually done in 
RMFT. With this prescription we were not able to reproduce the excitation 
energies of the just mentioned resonances. This has been possible only after 
introducing the repulsive Landau-Migdal (LM) delta force 

V LM (1,2) = g'l— J ti • r 2 (Ti ■ <r 2 5{n - r 2 ), (1) 
of the same magnitude (g f = 0.7) as the one used in the nonrelativistic 



calculations 1 1 



Here we wish to analyze whether the tensor (derivative) coupling of the p 
meson to the nucleon could generate a sufficiently repulsive nucleon-nucleon 
force in order to locate the GT resonance at the correct experimental energy 
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and in this way substitute the phenomenological LM force. The IA resonance 
is practically not affected by this part of the p-meson-exchange potential and 
therefore it will not be discussed so exhaustively as we do with the GT 
resonance. 

As mentioned above, it is not usual to include the tensor coupling of the 
vector mesons to the nucleon in RMFT. This is because its effect on the 
ground state is (rightly) thought to be small. On a more general perspective, 
however, there are two good reasons why one should do so. For one, according 
to the rules of effective field theory such terms should appear in the effective 
QHD Lagrangian [JT2 . For another, and perhaps more important reason for 



the phenomenological stand we are taking, it is well known that the tensor p- 
nucleon coupling gives a large contribution to the spin-isospin component of 



the nucleon-nucleon interaction ||13|| , and as such it could have an important 
effect on the dynamics of the GT resonance. 
Our Lagrangian density is now 



C = z/>(z 7/ ^-M)V> 

1 ... 1 o o I I 

2 M 2 - - -y- , 



1„ 1 



2 2 m^r 
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-R^ ■ R" u + -m p 2 p^ ■ p» - gpi/ry^Ti/) ■ p 



2M 



where 



- \f, v F^ - efaJ^p-W, (2) 



W = 9V - d^u", 

BT = SV-^-^/xp", (3) 
= d"A u -& / A lx . 



This Lagrangian is identical to that of reference | IT| , except for the p- nucleon 
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tensor coupling term (the one proportional to f p )P\ Therefore, following 
the same route one arrives at identical equations for the mean boson fields, 
except for that of the p meson (only the component p\ survives for spherical, 
definite-charge nuclei), which now takes the form 



(- 



-V 2 + m/) p% = g„p 3 (r) + ^ ■ PaM, (4) 



where the (vector) isovector density p% is as defined in [10| and we have 
introduced the tensor isovector density 

A 

p t3 = (^)icxT^ = J2UaiocT 3 U a . (5) 

a=l 

The summation is over all the occupied single-particle, positive-energy states 
U a , which obey the mean- field Dirac equation. This is also modified to 



{-ioc-V + (3[M + V s {r)]+V v {r) + {i^cx-r/r)V t {r)}U a = EJA a . 

(6) 

Again the scalar (V s ) and vector (V v ) potentials are as defined in [ID|, while 
the tensor potential, 

is the contribution from the tensor-coupling term in fl2|). 

The general structure and derivation of the RRPA for charge-exchange 
excitations, in the discretized spectral version we use, has been delineated in 
rpTOfl - An alternative, more detailed account can be found in [ 14| . The main 



ingredient is the residual interaction V. For a self-consistent calculation, 
this must be obtained from the same Lagrangian ([!]) used for the mean field. 
Also, since Fock terms are ignored in RMFT, we must consider only the direct 
matrix elements of V. Hence only the isovector mesons contribute, and we 
get V = Vk + V p , with, in the instantaneous approximation, 

K(l,2) = -(—) r 1 -r 2 {(T 1 -V 1 (T 2 -V 2 )Y{m^r 12 ), (8) 



lr There is a minor correction to be made in |10[ . One must replace g p by 2g p in eqs. 
(1) and (2) of that reference for consistency with the remaining equations. 
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VJ1,2) = Ti-T 2 



*-^°- V ) 1 («'-4« to - V ) 

-( S ,a + 4*xv).( S , £t+ A*xv) !j 



(9) 



where r 12 = \ri — r 2 | and F(m, r) = exp(— mr)/(47rr). 

For the numerical values of the parameters we follow mostly the philos- 
ophy of [fTC[| , adopting the parameter set NL1 [[| |T3j| . Yet, in view of the 
difficulties encountered by Ma et al. [[| in accounting for the El and EO 
giant resonances with the NL1 parameters, a few results for the TM1 model, 
worked out by Sugahara and Toki [[RJ, will be presented as well.Q Taking 
experimental values for the pion, the only new parameter is the p-nucleon 
tensor coupling constant f p . As mentioned in |T7]], the vector dominance 



model predicts for the ratio f p /g p = K p a value equal to the isovector mag- 
netic moment of the nucleon, i.e., fi p — \i n — 1 = 3.7. On the other hand, 
most meson-exchange models for the nuclear force use K p = 6.6 fI5 |. The 



former choice was preferred in the description of the ground state properties 



in closed shell nuclei within the relativistic Hartree-Fock approximation [[17 
Thus, the discussion that follows will mainly rely on the lower value for K p , 
even though we are aware of the fact that the inclusion of Fock terms can 
considerably change the adjusted values of the QHD parameters [[17J . 

Another point to consider is whether the inclusion of the tensor coupling 
term in the Lagrangian (0) does not sensitively affect the values of the re- 
maining parameters. Fortunately, while the contribution of this term is not 
strictly zero in RMFT, its effects on the single particle energies as well as on 
the ground state properties are certainly very small. We therefore feel justi- 
fied in keeping the remaining parameters fixed at their NL1 or TM1 values. 
With K p = 3.7, for instance, the spin-orbit splitting is modified in less than 
150 KeV. Similarly tiny effects on the energy per particle and the root- 
mean-square radii are displayed in Table [TJ. An interesting side remark can 

2 In the latter case, the w-meson self-interaction term, not appearing in (||), was also 
included in the numerical calculations. 

3 For identical particles, the NL1 paramerization yields significantly larger spin-orbit 
splittings than the TM1 model, while the opposite happens for nonidentical particles. 
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be made concerning the latter observables. It is well known that, at variance 
with the nonrelativistic calculations, it is a common feature of the relativistic 
models to overestimate the neutron skin thickness |T6|, |2B|. But, as seen from 
the results shown in Table [I], the tensor p-N coupling has the tendency to 
correct the RMFT for this handicap. This fact, in turn, could have very im- 
portant consequences on the estimates of the atomic parity nonconservation 

MM- 

The GT and IA resonances in 208 Pb were computed in RRPA, for both 
the NL1 and TM1 sets of parameters and within the same model space as 
that of [|Hj, i.e., including only OhQ and 2hQ excitations, and only those 
single-particle states that are bound at least for neutrons. For simplicity, 
we ignored the negative-energy states, although it has been shown || that 
they are required in principle even if the no-sea approximation is made for 
RMFT, since one needs a complete single-particle basis to develop a perfectly 
consistent RRPA. In fact, the transitions from Fermi- to Dirac-sea states are 
essential to ensure certain desirable features, such as current conservation 
and the removal of the spurious J n = 1~ translational state. However, such 
issues are not crucial for our present purposes and, furthermore, Ma et al. 



have shown in a recent calculation [26| that the contribution of the negative- 
energy states is of decisive importance only for the isoscalar modes. We 
therefore feel safe to leave their inclusion for a future, more sophisticated 
and detailed treatment of those isovector resonances. 

In Fig. |l] are shown the NL1 results for the GT strength distribution, 
both in terms of the individual strengths, 



E*2i(pIM 



\n) +J2 Y np(P\\ 



\n) 



pn 



np 



(10) 



and of a "strength function" obtained by replacing the spikes by Lorentzians 
of conveniently chosen widths A | IU|| , i.e., 



S(E) 



sx 



TT 



£a) 2 + A 2 



(11) 



where A"A and YA are, respectively, the forward and backward going RPA 
amplitudes for the state at excitation-energy E\. The upper, middle and 
lower panels correspond, respectively, to: (a) K p = 0,g' = 0; (b) K p = 
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3.7, g' = and (c) K p = 3.7, g' = 0.7. From these results one is induced to 
conclude that the tensor p-N coupling has a very small effect on the GT res- 
onance. That is, it seems as though this coupling could merely redistribute 
the GT strength in the energy region between 5 and 15 MeV, but was inca- 
pable of promoting it to the correct experimental energy. The latter is only 
achieved after introducing an LM force of the same magnitude that has been 
used in the previous calculation, where the just mentioned coupling has not 
been considered at all 
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The issue of the NN-force generated by the p-N 
coupling is, however, not so simple and deserves further discussion, which is 
presented below. Before proceeding, let us just mention that we have not 
noticed large differences between the NL1 and TM1 results for the IA and 
GT resonances. For instance, in the case (c) we get that these excitations 
are localized at: E IA (NL1) = 18.6 MeV and £ GT (NL1) = 19.5 MeV and 
Ei A (TMl) = 18.7 MeV and E GT (TM1) = 20.3 MeV, while the experimen- 
tal results are: E IA (exp) = 18.8 MeV and £ GT (exp) = 19.2 MeV. Thus, 
henceforth only the parametrization NL1 will be used. 

In the upper panel of Fig. |2] are confronted several diagonal J 71 = 1 + 
proton-particle neutron-hole matrix elements for the Vlm, V n , V p vv , V p VT 
and V p TT potentials. (The meaning of the upper indices is self-explanatory.) 
One can see, in particular, that the matrix elements of Vj T are very small 
in comparison with those coming from Vlm- However, when we rewrite Vj T 
in the form 



+ 




ti • T2P1P2 {-(a • V)i(a ■ V) 2 V(m p ,r 12 ) 



+ 



+ 1 Y(m p ,r 12 )S 12 



m 2 Y(m p , r 12 ) - 5(r 1 - r 2 ) 



(12) 



and evaluate different parts separately, we find out that the Yukawa and 
contact pieces in the last term engender, each one, very large matrix ele- 
ments. In fact, as shown in the lower panel, their individual values are larger 



than those of Vi 



LM; 



but the overall contribution to V p TT is small, 



because 



they cancel each other very strongly. (A similar cancelation, though not so 
pronounced, also occurs in the case of V n .) 

It should be remembered that the contact terms in V n and V p TT would be 
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smeared over a finite region if finite-nucleon-size effects (FNSE) were intro- 
duced, and they would be totally killed by realistic short range correlations 
(SRC). Yet, none of these two effects is considered in a mean field treatment, 
such as the present one. In return, it is common practice |5|, 17] to extract 



the contact parts from (J8|) and (H) by adding to the residual interaction the 
correction term 5V = 5V W + 5V P , with 

<$K(1,2) = U-) T X -T 2 <r x -* 2 5{r x -r 2 ), 



3 \m 7T/ 

SV P (1,2) = TrTjftft (a r a 2 + 2(r r ff 2 )i(ri-r 2 ) 



(13) 

For consistency, one must also perform such an extraction in the mean field 



part. Since, differently from the Hamiltonian formalism followed in [[17 
we are working within a Lagrangian formalism, we did this extraction in 
the baryon self-energy computed in the Hartree approximation (which is 
equivalent to RMFT). As a consequence the replacement V t — > V t + SV t has 
to be done in the Dirac equation (Q), with 




§Vt = l[AJL) 02 — _ r (14) 

being a correction that arises upon the extraction from the baryon self-energy 
of the contact part due to this derivative coupling in eq. (|j). But, when this 
recipe is implemented in the numerical calculation we get too much repulsion 
and the GT resonance is pushed up very high in energy. This comes from 
the fact that SV is basically a 5-force of the type (0), with 

= 1.6, (15) 

which is significantly larger than g' = 0.7. 

4 Note, however, that the contributions of the contact terms are nonzero when, both 
the FNSE, and the SRC are considered simultaneously pjj). 
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Note that in the nonrelativistic approximation the contact term also ap- 
pears in V p vv and V p VT , and instead of (p~5|) one would have 




\2MJ 



2.3. 



(16) 



There is no consensus on whether one should proceed in the same way in the 
relativistic case. Some authors exclude the contact terms only from V n and 
VT T [|T^| , while others do that for the full tt + p interaction 0, 28 1. That 



the potentials VY V and VY T also contain a contact term follows from the 



substitution [E9l 



2M 



(17) 



7m 

for the vector p-N coupling. 

It is worth noting that Toki and Weise |HJ have interpreted microscopi- 
cally the LM force as arising from the n + p meson-exchange model combined 
with the SRC and FNSE. In the static limit, which is used here, the result 
is IB1: 




g' LM (uJ = q = 0) = 

+ S^jr) [m) mg + mj' (18) 

where A is the cut-off mass for the pion-nucleon vertex and Wq 1 is the cor- 
relation length. For A = 1 GeV, m = m p and g p + f p = 17.2 this leads to 
g' LM (u = q = 0)= 0.67 0. (In the present work g p + f p = 23.4.) 



Our results can be summarized as follows: 

1. When the short range correlations are not considered, the tensor p- 
nucleon coupling plays only a minor role in the description of the GT 
resonances. 

2. If one tries to take these correlations into account by merely extracting 
the contact terms from the NN interaction, the GT resonance is pushed 
up too high in energy. 
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Thus, the simulation of the short range correlations by the simple-minded 
extraction of the contact terms alone is not a satisfactory procedure; at least 
not in the case of the heavier mesons. The explanation is that the contact 
terms in the 7r+p NN-interaction are not the only ones to be strongly modified 
by the short range correlations. In particular, because of the large p-meson 
mass, also the Yukawa terms generated in @ should be strongly reduced. 
We conclude hence that the implementation of, both realistic short range 
correlations, and finite-nucleon-size effects, in the context of the relativistic 
RPA, is required. Presently, we are working on this issue. 

Finally, let us mention that the tensor p-nucleon coupling plays an impor- 
tant role in the transverse spin response, and that some progress in assessing 
this through a relativistic many-body calculation has been made quite re- 
cently by Yoshida and Toki f32fl . 
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Figure 1: Gamow- Teller strength distribution for the parent nucleus 208 Pb for 
the parametrization NL1. The upper, middle and lower panels correspond, 
respectively, to: (a) K p = 0,g' = 0; (b) K p = 3.7, g' = and (c) K p = 
3.7, g' = 0.7. The spikes (r.h.s. scale) give the raw RRPA results and the 
continuous curve (l.h.s. scale), the strength function smoothed out by means 
of Lorentzians having widths of: (a) and (b) 3.0, and (c) 3.65 MeV. The 



strength function for the resonance peak extracted from experiment |22| is 
drawn in dotted line. 
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Figure 2: Diagonal matrix elements of: (a) the several terms of the 7r + p 
NN-interaction and (b) different pieces of V p TT , taken between proton-particle 
neutron-hole 1 + states in 208 Pb. The matrix elements of the Landau-Migdal 
contact force are also shown in both panels for comparison. The states are 
positioned at their unperturbed energies. 
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Table 1: RMFT results for the energy per particle and the root-mean-square 
radii of the neutron and proton point-particle distributions in 208 Pb com- 
puted with several values of the p-nucleon tensor coupling constant f p . The 
remaining parameters are kept fixed at their NL1 and TM1 values. 



K p = f p /g p 


E/A-M 


s](rl) 






[MeV] 


[fm] 


[fm] 


NL1 parameter set: 











-7.884 


5.795 


5.474 


3.7 


-7.882 


5.777 


5.480 


6.6 


-7.883 


5.763 


5.485 



TM1 parameter set: 






-7A 


374 


5.755 


5.485 


3.7 


-7A 


m 


5.741 


5.492 


6.6 


-7A 


m 


5.730 


5.497 


Experiment 


-7.8 


680 


5.593Q 


5.4520 


"Taken from Rcf. [ 
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b Taken from Ref. 
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c Taken from Refs. [£ 


| and H eq. (6.1)] 
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